--- layout: post title: Surface tension, capillary forces, and phase-volume relations date: Wed 03 Mar 2021 03:30:02 PM IST categories: ["Academic Writings", "Unsaturated Soil Mechanics"] --- -

1 Fundamental Concepts of Surface Tension

The forces between molecules of a fluid are called cohesive forces, the forces between the molecules of a fluid and a solid are called adhesive forces. The interaction of these forces results in surface tension and capillary, which are crucial for our discussion.

1.1 Surface tension

Figure 1: Surface tension as force per unit length along the surface

Through basic experiments of soap films on a wireframe, we can understand that surfaces exert force. We can further hypothesize that the liquid surface is in tension from this experiment, and we need to do work to increase the surface area.

A molecule in the midst of a liquid interacts with all its neighbor, i.e., it is attracted by the neighboring molecules equally in all the direction; However, a molecule on the surface has half of these interactions. Hence, the molecules on the surface are at higher energy states in comparison to internal molecules. The molecules of the liquid configure in such a way that the total energy is at minima. The surface energy of a liquid drop is proportional to number of surface molecules. Hence the liquid takes a configuration in which number of molecules on the surface are minimum, i.e., the surface area of the drop is minimum. To increase the surface area or to create new surface, additional energy will be required to bring the interior molecules on the surface. This energy will be equivalent to work required in breaking the bonds of interior molecules and bringing them on the surface where the bonds are half of the interior.

Suppose, we want to increase the surface area of a liquid by dAdA. The work required to bring the additional molecules on the surface will be proportional to dAdA (the additional number of molecules), i.e., we can write:

δW=γdA\delta{W}=\gamma{dA}

where γ\gamma is the surface tension - the energy required to increase the surface area by one unit.

The surface tension is also defined as force per unit length. Consider the soap film on the wire frame attached with a moving rod (fig. 1). Initially, the rod is in equilibrium; the rod is perturbed by small displacement, dxdx. The work done in this perturbation -

δW=Fdx2γLdx\delta{W}=Fdx-2\gamma{L}dx

Since the rod was in equilibrium -

δW=0\delta{W}=0 F=2γLdxF=2\gamma{L}dx

This suggests that the γ\gamma is the force exerted by the surface per unit length of the rod.

Typical values of surface tension tbl. 1 and a few remarks  [1]:

1.2 Laplace law

For a curved interface, the surface tension creates a pressure difference in fluids on both sides of the interface. In simple spherical drops, the pressure is higher on the concave side. This excess pressure has many consequences: formation of drops in a falling water jet, merging of smaller drops into larger ones, capillary adhesion (between two plates, between hairs, bewteen sand particles). In the following, we derive the pressure difference in terms of the surface tension.

1.2.1 Case 1 - spherical drop

Figure 2: Laplace pressure for oil drop in water

From fig. 2, the work required to increase the volume of drop

δW=po(4πR2δR)pw(4πR2δR)γow(8πRδR)\delta{W}=p_{o}\left(4{\pi}R^{2}\delta{R}\right)-p_{w}\left(4{\pi}R^{2}\delta{R}\right)-\gamma_{ow}\left(8{\pi}R\delta{R}\right)

At equilibrium, δW=0\delta{W}=0.

Δp=popw=2γowR\Delta{p}=p_{o}-p_{w}=\frac{2\gamma_{ow}}{R}

The above equation simply implies that for each spherical surface or interface, the pressure difference is 2γ12R\frac{2\gamma_{12}}{R}.

1.2.2 Case 2 - spherical bubble

Figure 3: Soap bubble - two interfaces or surfaces

A soap bubble is shown in fig. 3; there are two interfaces: inner and outer interface. For each interface, from inside to outside, the pressure increases by 2γsaR\frac{2\gamma_{sa}}{R}; hence

Δp=popi=4γsaR\Delta{p}=p_{o}-p_{i}=\frac{4\gamma_{sa}}{R}

1.2.3 Case 3 - generic surface

Figure 4: A part of generic surface

A small portion of a generic surface is shown in the fig. 4 with primary radii of curvature as R1R_{1} and R2R_{2}. In the following, we assume that the shape of surface is preserved for an arbitrarily small perturbation in which radius of curvature is increased by a small amount δz\delta{z}.

To preserve the shape, we need to ensure (from the fig. 4)

δx1R1=δx1+Δ(δx1)R1+δz\frac{\delta{x_{1}}}{R_{1}}=\frac{\delta{x_{1}}+\Delta{\left(\delta{x_{1}}\right)}}{R_{1}+\delta{z}}

similarly,

δx2R2=δx2+Δ(δx2)R2+δz\frac{\delta{x_{2}}}{R_{2}}=\frac{\delta{x_{2}}+\Delta{\left(\delta{x_{2}}\right)}}{R_{2}+\delta{z}}

The initial area is -

dA=δx1δx2dA=\delta{x_{1}}\delta{x_{2}}

Area after the perturbation -

dÃ=(δx1+Δ(δx1))(δx2+Δ(δx2))d\tilde{A}=\left(\delta{x_{1}}+\Delta{\left(\delta{x_{1}}\right)}\right)\left(\delta{x_{2}}+\Delta{\left(\delta{x_{2}}\right)}\right)

The increase in the area -

δdA=(1+δzR1)(1+δzR2)δx1δx2δx1δx2\delta{dA}=\left(1+\frac{\delta{z}}{R_{1}}\right)\left(1+\frac{\delta{z}}{R_{2}}\right)\delta{x_{1}}\delta{x_{2}}-\delta{x_{1}}\delta{x_{2}}

The work require to perturb the surface -

δW=piδx1δx2δz+poδx1δx2δz+γδdA\delta{W}=-p_{i}\delta{x_{1}}\delta{x_{2}}\delta{z}+p_{o}\delta{x_{1}}\delta{x_{2}}\delta{z}+\gamma\delta{dA}

At equilibrium -

δW=0\delta{W}=0

(pipo)δz=γ[δzR1+δzR2+(δz)2R1R2]\left(p_{i}-p_{o}\right)\delta{z}=\gamma\left[\frac{\delta{z}}{R_{1}}+\frac{\delta{z}}{R_{2}}+\frac{\left(\delta{z}\right)^{2}}{R_{1}R_{2}}\right]

Since δzR1R2\delta{z}\ll\ R_{1}R_{2},

(pipo)=γ(1R1+1R2)\left(p_{i}-p_{o}\right)=\gamma\left(\frac{1}{R_{1}}+\frac{1}{R_{2}}\right)

1.2.4 Capillary adhesion

Figure 5: Capillary adhesion between two plates

fig. 5 shows a drop of liquid mashed between two plates. The drop will take a form of funicular cylinder of radius RR and height HH. If we look from inside of the drop; there are two curvatures: concave one with radius R and convex one with radius H2cosθE\frac{H}{2\cos{\theta_{E}}}. Hence pressure difference between the liquid in the drop and outside gaseous medium is -

Δp=γ(1R1H2cosθE)\Delta{p}=\gamma\left(\frac{1}{R}-\frac{1}{\frac{H}{2\cos{\theta_{E}}}}\right)

As RR is sufficiently large -

Δp=2γcosθEH\Delta{p}=-\frac{2\gamma{\cos{\theta_{E}}}}{H}

The negative pressure acts as suction and attracts the plate towards each other by the force -

F=πR22γcosθEHF=\pi{R^{2}}\frac{2\gamma{\cos{\theta_{E}}}}{H}

Note - Rememeber that the Pascal law is only true when the surface is in equilibrium.

1.2.5 Minimal surfaces

1.2.5.1 Unstable water jet

A cylindrical jet of water of radius rr tend to break apart into drops; this can be explained using the Laplace law. As the cyliderical jet of water has radius rr; the pressure inside the drop is γr\frac{\gamma}{r} in access to atmospheric pressure. As there are no perfectly cylindrical jets, there are often irregularities which cause significantly high pressure and result in breakage of jet into drops. Further, the surface area of resulting drops should be less than the original surface area of cylindrical jet: the droplets form to minimze the surface area. Let us consider, the original jet of radius RR and length LL which resuls in nn droplets of radius rr. The ratio of original surface area to droplets surface area -

S0Sn=2πRL4πnr2\frac{S_{0}}{S_{n}}=\frac{2\pi{R}L}{4\pi{n}r^{2}}

From the conservation of volume (assuming density does not change) -

πR2L=43πnr3\pi{R^{2}}L=\frac{4}{3}\pi{n}r^{3}

This results in -

S0Sn=2r3R\frac{S_{0}}{S_{n}}=\frac{2r}{3R}

Since S0>SnS_{0}>{S_{n}} -

r>3R2r>\frac{3R}{2}

This suggests that the radius of droplets will be higher than the original jet.

1.2.5.2 Drop on a fibre

Figure 6: Minimal surface - a liquid drop on the fibre

Consider a drop of liquid lying on a fibre pipe of radius bb; the contact angle at the interface is zero, i.e., the liquid is able to wet the fibre. From the Pascal law at a point on the surface of the drop (refer fig. 6) -

Δp=γ(1zcosθz̈(1+ż2)3)\Delta{p}=\gamma\left(\frac{1}{\frac{z}{\cos{\theta}}}-\frac{\ddot{z}}{\sqrt{\left(1+\dot{z}^{2}\right)^{3}}}\right)

The equation of equilibrium at each point of the surface of drop is -

(1z1+ż2z̈(1+ż2)3)=Δpγ\left(\frac{1}{z\sqrt{1+\dot{z}^{2}}}-\frac{\ddot{z}}{\sqrt{\left(1+\dot{z}^{2}\right)^{3}}}\right)=\frac{\Delta{p}}{\gamma}

Further, we can derive a first order ordinary differential equation representing the surface of the drop by considering the equilibrium of a part of the drop, as shown in fig. 6.

2πγb+2πγzcosθπ(z2b2)Δp=0-2\pi\gamma{b}+2\pi\gamma{z}\cos{\theta}-\pi\left(z^{2}-b^{2}\right)\Delta{p}=0

zcosθb(z2b2)Δp2γ=0z\cos{\theta}-b-\left(z^{2}-b^{2}\right)\frac{\Delta{p}}{2\gamma}=0

z1+ż2b(z2b2)Δp2γ=0\frac{z}{\sqrt{1+\dot{z}^{2}}}-b-\left(z^{2}-b^{2}\right)\frac{\Delta{p}}{2\gamma}=0

At z=Lz=L, ż=0\dot{z}=0 -

Lb(L2b2)Δp2γ=0L-b-\left(L^{2}-b^{2}\right)\frac{\Delta{p}}{2\gamma}=0

The excess pressure -

Δp=2γ1L+b\Delta{p}=2\gamma\frac{1}{L+b}

1.3 Wetting - contact between three phases

Figure 7: Wetting phenomena

Wetting is the phenomena of spreading of liquid on a solid substrate. Wetting can be characterized into -

  1. Total wetting - liquid has strong affinity to solid substrate.

  2. Partial wetting - otherwise.

These caan be identified using the spreading parameter.

1.3.1 Spreading parameter

The spreading parameter is difference in the surface energy per unit area of interface before and after the deposition of the drop.

S=γSG(γSL+γLG)S=\gamma_{SG}-\left(\gamma_{SL}+\gamma_{LG}\right)

Total wetting occurs when the surface energy decreases, i.e., S>0S>0; else partial wetting happens with S<0S<0. In total wetting, the drop completely spreads over the surface in other to minimize the surface energy; whereas, in partial wetting, the drop forms a curved cap at equilibrium resulting in a contact angle θE\theta_{E}. We can further say -

  1. θE=0\theta_{E}=0 represents total wetting.
  2. θE<π2\theta_{E}<\frac{\pi}{2} \rightarrow partial wetting \rightarrow mostly wetting
  3. θE>π2\theta_{E}>\frac{\pi}{2} \rightarrow partial wetting \rightarrow mostly non-wetting

1.3.2 Law of Young-Dupre

Figure 8: Determination of contact angle

Using the force balance at line of contact or triple line, we can obtain a relation for contact angle θE\theta_{E} as

γSG(γSL+γLGcosθE)=0\gamma_{SG}-\left(\gamma_{SL}+\gamma_{LG}\cos{\theta_{E}}\right)=0

γLGcosθE=γSGγSL\gamma_{LG}\cos{\theta_{E}}=\gamma_{SG}-\gamma_{SL}

We can further write -

S=γLG(cosθE1)S=\gamma_{LG}\left(\cos{\theta_{E}}-1\right)

or

cosθE=1+SγLG\cos{\theta_{E}}=1+\frac{S}{\gamma_{LG}}

Clearly, θE\theta_{E} is not defined when S>0S>0.

We can also derive this relation by using stationarity condition for WW at equilibrium.

δW=(γSGγSL)dxγLGdxcosθE\delta{W}=\left(\gamma_{SG}-\gamma_{SL}\right)dx-\gamma_{LG}dx\cos{\theta_{E}}

At equilibrium -

δW=0\delta{W}=0

1.4 Capillary

The capillary action or capillarity is the phenomena of a liquid flowing into narrow cavities or channels against the direction of gravity due to interfacial tension. Read the excerpt from  [2] about the initial investigation on capoillarity -

with a view to explain arise apparent anomalies. After long groping in the dark, it was found to be desirable to discover by experiment what were the actual phenomena which required explanation. Hawksbee found that the height to which a fluid would rise in From a capillary tube of given radius was the same for all . thicknesses of the tube. was apparent that the attracting force of the tube was situated at or near the inner surface of the tube. But he does not appear to have taken account of the mutual attractions of the particles of the fluid. Jurin also found that the height of the column of fluid supported by capillary action depended solely upon the interior diameter of the tube at the upper surface of the fluid. From this he con- cluded that the column of fluid raised’ by Capillary Action was supported by the attraction of the periphery or section of the tube to which the upper surface of the fluid cohered or was contiguous.

Clairaut was the first to attempt to explain capillary phenomena on right principles, by referring them to the mutual attraction of the particles of the fluid, and to the attraction of the particles of the solid on the particles of the fluid; and supposing these attractions to depend upon the same function of the distance, he concludes that even if the attraction of the capillary tube be of a less intensity than that of the water, provided the intensity of the latter attraction be not twice as great as that of the former, the water will still rise in the tube (p. 121). Clairaut supposed that the attraction was sensible only at very small distances (p. 113).

Laplace (1706) realized that the liquid rise in a capillary tube if the dry surface energy of the tube, γSG\gamma_{SG}, is greater than wet surface energy, γSL\gamma_{SL}, per unit area  [3]. Similar to the spreading parameter, an imbibition parameter, I, is defined as

I=γSGγSLI=\gamma_{SG}-\gamma_{SL}

I=γLGcosθI=\gamma_{LG}\cos{\theta}

when I is positive \rightarrow liquid rises, otherwise it drops in the tube.

1.4.1 Capillary length

Near a miniscus, the capillary force dominates till the capillary length after which gravity dominates. Capillary length can be estimated by -

ρgκ1=γκ1\rho{g}\kappa^{-1}=\frac{\gamma}{\kappa^{-1}}

κ1=γρg\kappa^{-1}=\sqrt{\frac{\gamma}{\rho{g}}}

Figure 9: Capillary length

From the fig. 1, we can derive the following -

pApatm=γz̈(1+ż2)3p_{A}-p_{atm}=-\gamma\frac{\ddot{z}}{\sqrt{\left(1+\dot{z}^{2}\right)^{3}}}

At the capillary length, ż1\dot{z}\ll{1} -

γz̈ρgz=0\gamma\ddot{z}-\rho{g}z=0

z=z0exp(κx)z=z_{0}\exp{\left(-\kappa{x}\right)}

1.4.2 Capillary rise in a tube

We can find the capillary rise using the Laplace pressure law. The pressure near the meniscus on the liquid side of the interface is -

p=patm2γcosθERp=p_{atm}-\frac{2\gamma\cos{\theta_{E}}}{R}

also

p=patmρgHp=p_{atm}-\rho{g}H

From above equations -

H=2γRcosθEρgH=\frac{2\gamma{R}\cos{\theta_{E}}}{\rho{g}}

This can also be derived using equilinrium of forces and by minimization of energy.

Figure 10: Capillary rise in a tube

1.5 Capillary between two particles

Figure 11: Capillary bridge between spherical particles
Figure 12: Equilibrium of a section of capillary bridge

fig. 11 shows a capillary bridge between two spherical particles. Understanding of capillary between two particles is important in understanding the flow behaviour in unsaturated medium, tensile strength, and moisture retention of unsaturated soils. The shape of capillary bridge is unknown apriori and is a function of contact angle, separation between the spheres (representing the roughness), surface tension, and volume of the liquid. This problem can be solved using the minimization of surface area with the given volume and further constraints on shape of sphere and contact angle. Here we will attempt the solution using the equilibrium method (we could also solve it using the Laplace law).

Consider a section ABCD of the capillary bridge shown in fig. 5. As the capillary bridge in equilibrium any part of it will also be in equilibrium -

[2πy0]γ+[2πy]γcosα+[πy2πy02]ψ=0-\left[2\pi{y_{0}}\right]\gamma+\left[2\pi{y}\right]\gamma\cos{\alpha}+\left[\pi{y^{2}}-\pi{y_{0}}^{2}\right]\psi=0

Substituting cosα=11+ẏ2\cos{\alpha}=\frac{1}{\sqrt{1+\dot{y}^{2}}} and simplifying further, we get -

dydx=y2(y0+ψ2γ(y02y2))21\frac{dy}{dx}=\sqrt{\frac{y^{2}}{\left(y_{0}+\frac{\psi}{2\gamma}\left(y_{0}^{2}-y^{2}\right)\right)^{2}}-1}

We can solve the above ordinary differential equation numerically using the boundary condition  [4]:

  1. x=0x=0, y=y0y=y_{0}, dydx=0\frac{dy}{dx}=0.
  2. x=s+R(1cosβ)x=s+R\left(1-\cos{\beta}\right), y=Rsinβy=R\sin{\beta}, dydx=cot(β+θ)\frac{dy}{dx}=\cot{\left(\beta+\theta\right)}.

The second boundary condition can be used to determine the capillary sunction ψ\psi.

1.6 Application of surface tension

  1. Growth of plants  [5]
  2. Sandcastles  [6]

2 Phase properties and relations

2.1 Phase properties

2.1.1 Density and specific volume

Density of the soil particles is given as

ρs=MsVs\rho_{s}=\frac{M_{s}}{V_{s}}

Density of the air phase is given by

ρa=MaVa\rho_{a}=\frac{M_{a}}{V_{a}}

The specific gravity is defined as the ration of the density of the solid and density of air at 44^\circC and at 1 atm pressure

Gs=ρsρwG_{s}=\frac{\rho_{s}}{\rho_{w}}

Air behaves as mixture of several gases and also varying amounts of water vapour.

Air is found in two states - Dry air (without water vapour) - Moist air (with water vapour)

Air behaves like an ideal gas and the behaviour is governed by

uaVa=MawaRT\bar{u}_{a}V_{a}=\frac{M_{a}}{w_{a}}RT

where, ua\bar{u}_{a} is the absolute air pressure, VaV_{a} is the volume of the air, MaM_{a} is the mass of the air, waw_{a} is the molecular mass of the air, RR is the universal gas constant, TT is the absolute temperature.

In an isolated system, air follows the Boyle’s law

ua1Va1=ua2Va2\bar{u}_{a_{1}}V_{a_{1}}=\bar{u}_{a_{2}}V_{a_{2}}

Density of the air is given by

ρa=MaVa=uawaRT\rho_{a}=\frac{M_{a}}{V_{a}}=\frac{\bar{u}_{a}w_{a}}{RT}

The concentration of the water vapour in the air is commonly expressed in terms of the relative humidity

RH=100uvuv0RH=\frac{100\bar{u}_{v}}{\bar{u}_{v0}}

where, ua\bar{u}_{a} partial pressure of water vapour at the same temperature, uv0\bar{u}_{v0} is the saturation pressure of water vapour at the same temperature.

The density of air decreases as the relative humidity increases i.e. moist air is lighter than the dry air.

2.1.2 Viscosity

2.2 Interaction of air and water

2.2.1 Solid, liquid, and vapour states of water

2.2.2 Water vapour

2.2.3 Air dissolving in water

2.3 Volume mass relations

2.3.1 Porosity

2.3.2 Void ratio

2.3.3 Degree of saturation

2.3.4 Water content

2.3.5 Soil density

2.3.6 Basic mass volume relationship

2.3.7 Changes in mass volume properties

2.3.8 Density of mixture subjected to compression

3 Readings

  1. De Gennes, P.G., Brochard-Wyart, F. and Quéré, D., 2013. Capillarity and wetting phenomena: drops, bubbles, pearls, waves. Springer Science & Business Media.

  2. Chemistry: molecular approach - liquids, solids, and molecular forces

4 References

[1]
L. Trefethen, Surface Tension in Fluid Mechanics, Lubricating Oil 25, 35 (1969).
[2]
F. Bashforth and J. C. Adams, An Attempt to Test the Theories of Capillary Action by Comparing the Theoretical and Measured Forms of Drops of Fluid (University Press, 1883).
[3]
P.-G. De Gennes, F. Brochard-Wyart, and D. Quéré, Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves (Springer Science & Business Media, 2004).
[4]
F. Molenkamp and A. H. Nazemi, Interactions Between Two Rough Spheres, Water Bridge and Water Vapour, Géotechnique 53, 255 (2003).
[5]
I. Woodward, Tall Storeys, Nature 428, 807 (2004).
[6]
D. J. Hornbaker, R. Albert, I. Albert, A.-L. Barabási, and P. Schiffer, What Keeps Sandcastles Standing?, Nature 387, 765 (1997).